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ABSTRACT 

In this paper, in the framework of the secondary infall model, the correlation 
between the central surface density and the halo core radius of galaxy, and cluster 
of galaxies, dark matter haloes was analyzed, this having recently been studied on a 
wide range of scales. We used Del Popolo (2009) secondary infall model taking into 
account ordered and random angular momentum, dynamical friction, and dark matter 
(DM) adiabatic contraction to calculate the density profile of haloes, and then these 
profiles are used to determine the surface density of DM haloes. The main result is 
that (the halo characteristic radius) is not an universal quantity as claimed by 
Donato et al. (2009) and Gentile et al. (2009). On the contrary, we find a correlation 
with the halo mass M200 in agreement with Cardone & Tortora (2010), Boyarsky 
at al. (2009) and Napolitano et al. (2010), but with a significantly smaller scatter, 
namely 0.16 ± 0.05. We also consider the baryon column density finding this latter 
being indeed a constant for low mass systems such as dwarfs, but correlating with 
mass with a slope a — 0.18±0.05. In the case of the surface density of dark matter 
for a system composed only of dark matter, as in dissipationless simulations, we get 
a — 0.20 ± 0.05. These results leave little room for the recently claimed universality 
of (dark and stellar) column density. 
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1 INTRODUCTION 

The ACDM model is remarkably successful in fitting a wide 
range of data on cosmological scales, from the Hubble dia- 
gram of Type la Supernovae (Kowalsky et al. 2008) to the 
matter power spectrum with its Baryonic Acoustic Oscilla- 
tion features (Percival et al. 2010) and the anisotropy and 
polarization spectrum of cosmic microwave background ra- 
diation (Komatsu et al. 2011). Notwithstanding this long 
list of successes awarding the ACDM scenario the name 
of concordance model, there is still some tension at scales 
from a few kpc to tens of pc so that great amount work 
has been recently dedicated to test what ACDM predicts 
for galaxies and clusters. On the one hand, numerical CDM 
simulations typically form a number of satellite haloes, i.e. 
small haloes and galaxies orbiting around a massive central 
system, which is far larger than what is indeed observed. 
As a clear evidence of this substructure problem (Moo et 
al. 1999), one can note that the best studied satellite sys- 
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terns of the Milky Way and Andromeda galaxies are hard 
to reconcile with the predicted abundance of small -mass 
dark matter clumps even if one takes into account detec- 
tion effects. On the other hand, rotation curves of spiral 
galaxies favor cored dark halo models in contrast with the 
cusps predicted by the ACDM model de Blok et al. (2009), 
while a severe discrepancy is also observed for what con- 
cerns the size and the angular momenta of galactic disks 
(Mayer et al. 2008). Although both theoretical models and 
increasingly accurate observations have been put down dur- 
ing the past 10-15 years, the quoted problems still stay on 
the ground. Finding the origin of such discrepancies is com- 
plicated mainly due to the lack of a well understood method 
to describe the complex hydrodynamical phenomena in hot 
dense plasma and physical interactions of baryons and ra- 
diation which follow the initial (well understood) growth 
and collapse of the initial density perturbations. As a valu- 
able help to address this problem, one can look for scaling 
relations among DM haloes parameters and stellar quanti- 
ties in order to better constrain the formation scenarios and 
the DM properties. In this context, an intriguing property 
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of dark matter haloes was noted by Kormendy & Freeman 
(2004) based on halo parameters obtained by mass model- 
ing 55 spiral galaxy rotation curves within the framework 
of the Maximum Disk Hypothesis (MDH) . Among other re- 
lations between the halo parameters, they found that the 
quantity (iod = pofcQ, proportional to the halo central sur- 
face density for any cored halo distributions, is nearly inde- 
pendent of the galaxy blue magnitude. In particular, they 
found that this quantity takes a value of ~ IOOMq /pc~ 2 . 
Donato et al. (2009, D09), by means of rotation curves of 
~ 1000 spiral galaxies, the mass models of individual dwarf 
and spiral galaxies and the weak lensing signal of elliptical 
and spirals, found strong evidences for the constancy of the 
central DM column density, over 12 orders of magnitude in 
luminosity. They found log /ion = 2.15 ±0.2, in units of 
log(M©/pc 2 ). Gentile et al. (2009, G09) extended the result 
to the luminous matter surface density, finding that total 
luminous-to-dark matter ratio within one halo scale-length 
is constant. Opposite results were obtained by Napolitano, 
Romanowsky & Tortora (2010, NRT10) that showed that, 
on average, the projected density of local early type galax- 
ies (ETGs), within effective radius is systematically higher 
than the same quantity for spiral and dwarf galaxies, point- 
ing to a systematic increase with halo mass as suggested by 
Boyarsky et al. (2009a, B09). These last extended the sam- 
ples above analyzed to both group and cluster scale systems 
finding that the dark matter column density, S, (defined in 
Sect. 3 and equivalent to /j,qd in the case of fit with Burkert 
profile (Burkert 1995)) is given by 

lo g 5 = 0.21 log ^|^ + 1.79 (1) 

with S in Mq/pc 2 . Cardone & Tortora (2010) (CT10), by 
modeling the dark halo with a Navarro - Frenk -White pro- 
file and assuming a Salpeter initial mass function (IMF) to 
estimate stellar masses, found that the column density and 
the Newtonian acceleration within the halo characteristic 
radius r s and effective radius R e f / are not universal quanti- 
ties, but correlates with the luminosity Lv , the stellar mass 
M, and the halo mass -M200, contrarily to what is expected 
from D09 and G09, and in agreement with B09. 

In order to try to discriminate among these results and 
find an explanation and analytical derivation of the surface 
density of halos, we analyze the problem using the secondary 
infall model (SIM) introduced in Del Popolo (2009) (here- 
after DP09), taking into account ordered and random an- 
gular momentum, dynamical friction, and baryon adiabatic 
contraction. The plan of the paper is the following. In Sect. 
2, we summarize the model. In sect. 3, we discuss the results 
and Sect. 4 is devoted to conclusions. 



2 MODEL 

The above discussed results about a constant surface den- 
sity (D09; G09) of DM together with the ones suggesting a 
mass dependence (B09; CT10) are substantially grounded on 
fitting observational properties. With the exception of Bo- 
yarsky et al. (2010) (B10), no qualitative explanation and/or 

po and ro are, respectively, the central density and core radius of 
the adopted pseudo-isothermal cored dark matter density profile. 



analytical derivation of the quoted results has been proposed 
so far. B10 used a SIM and supposed that the halo density 
profile evolves in a self-similar way as previously done by 
Fillmore & Goldreich (1984) (FG84), Bertschinger (1985) 
(B85) and Skivie et al. (1997) (S97). They arrived at the 
conclusion that S oc M 1 ^ 3 t~ i ^ 3 qualitatively in agreement 
with B09 simulations and CT10 results, in the sense that S 
increases with mass, but the slope is much bigger than what 
obtained by B09 and CT10. 

The quoted discrepancy could be connected to the over- 
simplified structure of the SIM used. Even if they tried to 
use the improved SIM of S97, taking account angular mo- 
mentum, these effects on collapse, in reality, are not properly 
taken into account. 

In the following, we will shortly describe how imple- 
menting the SIM approach in oirder to obtain the dark mat- 
ter halo density profile. The interested reader can find more 
in DP09. 

The halo profiles are derived by using the improved sec- 
ondary infall model of DP09. In the model, a bound mass 
shell of initial comoving radius Xi expands up to a maxi- 
mum radius (or turnaround radius) x ta ,. As successive shells 
expand, they acquire angular momentum and then contract 
on orbits determined by the angular momentum itself, while 
dissipative processes and eventual violent relaxation inter- 
vene to virialize the system converting kinetic energy into 
random motions. The final density profile may then be com- 
puted as : 



p(x) = 



(x/x t!i ) 



1 + 



din (x/x t!l 
din a; ta 



(2) 



with pta.{xta) the density at turnaround and x/x ta , referred 
to as the collapse factor (see Eq. A18 in DP09). To describe 
the proto-haloes density profile, DP09 considered the profile 
of a peak in the density field generated according to the 
matter power spectrum of Bardeen et al. (1986) and then 
took into account angular momentum, dynamical friction 
and the presence of baryons as described below. 

First, the angular momentum is decomposed in an or- 
dered component, related to the tidal torques experienced 
by proto-haloes, and a random component connected to ran- 
dom velocities (Ryden & Gunn 1987). The ordered term is 
computed following Ryden ( 1988a, b), while the random part 
is assigned to proto - structures according to Avila - Reese et 
al. (1998) and Ascasibar et al. (2004). A term related to the 
dynamical friction force has been explicitly introduced in 
the equations of motion and evaluated Kandrup (1980) di- 
viding the gravitational force into an average and a random 
component generated by the clumps in the hierarchical uni- 
verse. Finally, adiabatic contraction of the halo, due to the 
baryonic collapse, has been taken into account through the 
formalism of Klypin et al. (2002) and Gnedin et al. (2004), 
also including the exchange of angular momentum among 
baryons and dark matter (Klypin et al. 2002, and Gnedin et 
al. 2004). 

As previously reported, the details of the model are 
given in DP09. The profile formation is due to a complex 
baryon-dark matter (DM) interplay that can be summarized 
as follows 

Initially the proto-structure is in the linear phase, it ex- 
pands, reach a maximum of expansion and then collapse. 
Baryons trapped inside the potential wells of DM halos are 
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subject to radiative dissipation processes which give rise to 
clumps and self-gravitating clouds before it collapses to the 
halo center and condenses into stars and galaxies. The stage 
of baryons cooling and stars formation happens as described 
in Ryden (1988) (Sect. 4: baryionic dissipation). In the in- 
fall, the baryons compress the DM halos (AC), producing a 
steepening of the DM density profile. When clumps reaches 
the central high density regions they experience a dynami- 
cal friction force from the less massive DM particles as they 
move through the halo. Dynamical friction acts as an an- 
gular momentum engine (Tonini, Lapi & Salucci 2006), and 
energy and angular momentum are transferred to DM (El- 
Zant et al. 2001), increasing its random motion, and giving 
rise to a motion of DM particles outwards. Moreover, or- 
dered angular momentum mainly acquired in the expansion 
phase, through tidal torques, gives rise to nonradial motions 
in the collapse phase that amplifies the effects of the previ- 
ous mechanism. Then the joint effect of angular momenta 
(ordered and random) and dynamical friction trasporting 
it from baryons to DM, overcomes that of the AC and the 
profile starts to flatten. 

The comparison of the quoted model with SPH simula- 
tions has been performed in several papers (e.g., Del Popolo 
2012, Fig. 3), showing a good agreement with hydrodinam- 
ical simulations. 

The final product of this halo formation method gives 
the DM density profile as function of the radius r and the to- 
tal halo mass M v i r - The latter is the only parameter needed 
in order to specify the halo density, being the halo inner 
slope a function of the virial mass as well. Such a depen- 
dence breaks the universality of the halo profiles and favor 
Burkert models at dwarf scales, and models steeper than 
NFW at normal galaxy scales. 



3 RESULTS 

3.1 The S-M relation 

D09 and G09 analyzed rotation curves and weak lensing data 
for a sample of dwarf, spiral and elliptical galaxies fitted by 
the Burkert profile: 

= i X TL ^ ( 3 ) 
[r + r )(H + rg) 

The fit to the rotations curves yields the values of the two 
structural DM parameter (i.e., ro and po), then the surface 
density is calculated as fioD = po r o- The previous approach 
is based on the assumption that the galaxies studied, namely 
dwarf, spiral and elliptical galaxies can be fitted by cored 
models Q Now, the D09 and G09 samples contains dwarfs, 
spirals and ellipticals. 

For what concerns dwarfs, a large part of the studies of 
their density profiles indicate that the inner part of density 



profiles is characterized by a core- like structure (e.g., Flores 
& Primak 1994; Moore 1994; Kravtsov et al. 1998; de Blok 
& Bosma 2002; de Blok, Bosma & McGaugh 2003; Gentile 
et al. 2004, 2007; Blaise-Ouelette et al. 2004, Spano et al. 
2008, Kuzio de Naray et al. 2008, 2009, and Oh et al. 2010), 
but some studies found that density profiles are cuspy. Other 
studies concluded that density profiles are compatible with 
both cuspy and cored profiles (van den Bosch et al. 2000; 
Simon et al. 2005; de Blok et al. 2008). Very interesting are 
the results of Simon et al. (2005), and de Blok et al. (2008). 
The first showed that the power-law index, a, describing 
the central density profile, span the range from a — (NGC 
2976) to a = 1.28 (NGC 5963), with a mean value a = 0.73. 
Del Popolo (2012), showed how cuspy profiles can form even 
in the case of dwarfs. The second (de Blok et al. 2008), by 
using the THINGS sample, found that for low mass galaxies 
a core dominated halo is clearly preferred over a cusp-like 
halo, while for massive, disk dominated galaxies, the two 
models fit apparently equally well. In other terms, while a 
large part of dwarfs are well fitted by core-like profiles, some 
of them are not (see also Del Popolo 2012). If this cannot be 
taken as an evidence in favour of cusps, it shows that using 
core profile for every spiral is not always correct, and if one is 
using a sample containing also ellipticals the problem is even 
more serious. Moreover, THINGS shows that there is a trend 
to more cuspy profiles going to higher masses. The trend is 
confirmed by the fact that ellipticals are characterized by 
cuspy profiles and similarly clusters of galaxies have mainly 
cuspy profiles, (see the following of the discussion). 
NRT10 studied the DM surface density finding that their 
result is in good agreement with those of D09, namely they 
found a constant surface density on average in the mass 
range going from dwarfs to late types, but differently from 
D09 they found a systematic increase with mass as found 
by Boyarsky et al. (2009). Moreover, according to them, the 
consistency of the NFW models with the dwarf and spiral 
data shown in their Fig. 9 suggests that the surface density 
result does not necessarily imply cored DM haloes. 

In the case of ellipticals, the situations is even more 
complicated. Several studies (e.g., Mamon & Lokas 2005) 
showed that the DM profile of ellipticals is cuspy. NRT10 
found that early-type-galaxie (ETG) violate the constant 
density scenario for the other galaxies by a factor of ~ 10 
on average, and a factor of ~ 5 in the same mass regime, in 
agreement with Gerhard et al. (2001), Thomas et al. (2009), 
and Boyarsky et al. (2009). 

The quoted discussion imply that when we study the 
surface density it is preferable to have a "method" that is 
able to fit not only cored density profiles but also cuspy or 
intermediate among the two types. 

Alternatively we can determine the surface density, in 
a more general way than that used by D09 and G09, by 
means of a dark matter column density, which represents 
an average over the central part of the object: 



2 The existence of a constant central surface density of DM in 
galaxies is independent from the core model used (D09). The 
choice of of a Burkert profile (which is a model with core) is 
due to the fact that at radii > 0.3r v i r (where r v i r is the virial 
radius) the Burkert profile and the NFW profile converge, while 
at small radii and for appropriate values of ro it reproduce the 
NFW profile. 



S = ~S J rdr J dz PDM{sJr 2 + z 2 ) (4) 

The integral over z is extended to the vitial of the halo. 
Definition Q has as a consequence that S scales as the 
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dark matter surface density within r* (S oc p*r*), but 5 
is more universal, being denned for any density profile (not 
necessarily in the core). In other terms, if the profile is cored 
using S we will obtain the same result for the surface density 
of the Burkert's fit, while in the case the profile is not cored 
S will give a more precise value for the surface density than 
using the Burkert's model fit. 

The quoted method was already used by B09 who ex- 
tended the analysis of D09 and G09 to galaxies and galaxy 
clusters and fitted the DM profiles by means of three DM 
profiles models, namely Burkert profile (Eq. [3|, pseudo- 
isothermal (ISO) profile 



p(r)iso = 



1 +r 2 /V2 



(5) 



(6) 



where r c is the ISO core radius, and NFW profile: 

P(r)NFW = -7— ^2 

r(l + r/r s ) 

where r s is the typical NFW characteristic radius. 

CT10 used two DM profiles for the fit, namely Burkert 
profile and ISO. 

In order to determine S, we use Eq. Q together with 
the density profiles obtained with the model of the previ- 
ous section, which are plotted in Fig. 1, 3, and 5 of DP09. 
As there shown, density profiles become monotonically shal- 
lower inwards, down to the innermost resolved point, with no 
indication that they approach power-law behavior. This is in 
agreement with recent N-body simulations (e.g., Navarro et 
al. 2010). Navarro et al. (2010) showed that density profiles 
deviate slightly but systematically from the NFW model, 
and are approximated more accurately by a fitting formula 
where the logarithmic slope is a power-law of radius, the 
Einasto profile 



ln(p(r)/p_ a ) = -2/a[(r/r_ a ) a - 1], 



(7) 



where r_2 and p_2 are connected to the scaling radius and 
density of the NFW profile (see the following) by r_2 = 
r s and p_2 = Pa/4. A comparison between the results of 
our model and dissipationless N-body simulations and SPH 
simulations are shown in Del Popolo (2011). 

In Fig. 1, S(r s ) vs M200 is plotted for our model and for 
different authors'. From bottom to top: the best fit linear re- 
lation, found using the D' Agostini (2005) fit, is represented 
by the dotted line, while the solid line is related to the di- 
rect fit methods obtained in CT10. The long-short dashed 
line represents S(r a ) found in the framework of the model 
of this paper. The short-dashed line is connected to the 
prediction and results of the ACDM N-body simulation by 
Maccio et al. (2008). The long-dashed line is related to the 
B09 best fit linear relation. The dot-dashed line is related 
to the secondary-infall model (B10) theoretical prediction. 

The maximum likelihood fit for the correlation 
log 5(r s )-log A/200 in this model is 



3 Parameters of different profiles that fit the same DM density 
distribution are related (for example, r s for NFW is equal to 6.1r c 
for ISO and equals to 1.6t\b for BURK). Choosing these values as 
r* in each case, one finds that the value of S for NFW and ISO 
differ by less than 10% (the difference between NFW and BURK 
is ~ 2%), see Supporting Information for details. 
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Figure 1. S(r a ) as a function halo mass M20O' From bottom to 
top we have that: the dotted line represents the best fit linear 
relation using the D' Agostini (2005) fit, while the solid line the 
direct fit methods, obtained in CT10. The long-short dashed line 
represents S(r s ) obtained with the model described in the present 
paper. The short-dashed line represents prediction and the results 
from the ACDM N-body simulation of Maccio et al. (2008). The 
long-dashed line represents the B09 best fit linear relation. The 
dot-dashed line represents the secondary-infall model (B10) pre- 
diction. Black points and open boxes plot S(r s ) as a function halo 
mass M 2 oo, respectively for the NFW+Salpeter model and when 
a Chabrier IMF is used, obtained by CT10. 



log S(r s ) = 0.16 log( 



M200 



' 10 12 M Q 
while in the case of CT10 is 

log<S(r s ) =0.161og( 



+ 2.23 



+ 2.11 



(8) 



10 12 M Q 

and a steeper slope is obtained with the stellar mass, M. 
M 200 



(9) 



log S{r s ) =0.291og( 



10 n M ( 



+ 2.17 







(10) 



The marginal constraints on the scaling relationships pa- 
rameters for the correlation log <S(r s )-log M200 is 0.161q'q5 
in this model, while for CT10, considering a fiducial 
NFW+Salpeter model, is O.M+Qjg and for the correlation 
logS(r s ) - logM, is 0.29±°;^ (Table 3 of CT10). 

The best fit concerning DM halos in CT10 is shallower 
than the B09 and the result of the present paper, altough 
the slope is consistent with B09 (0.21) and the present pa- 
per within the large error bars. One has to recall that while 
CT10 have explicitly taken into account the correlated er- 
rors on S(r a ) and M 2 oo, it is not known if B09 does the same. 
Then, the difference in slope is only the outcome of the use 
of different algorithms on noisy data. However, if we com- 
pare the results of the present paper, those of B09 (over the 
same mass range) the estimates from the ACDM N- body 
simulations of isolated haloes from Maccio, Dutton & van 
den Bosch (2008) and the predictions from the secondary in- 
fall model (B09), are on average larger than those of CT10 
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Figure 2. S(r s ) as a function halo mass A/200- From bottom to 
top we have that: the short-dashed line represents the prediction 
and the results from the ACDM N-body simulation of Maccio et 
al. (2008); the long— dashed line the B09 best fit linear relation; the 
solid line represents the direct fit methods, obtained in CT10. The 
short-long dashed line represents S(r s ) obtained with the model 
described in the present paper. 



S(rs) values. A possible reason for the difference that the 
larger values in literature neglecting the stellar component. 

Another fundamental point to remark is the fact that in 
CT10, because of the large errors the log <S(r s )-log A/200 is 
consistent even with a slope 0. In our paper, because of the 
smaller errors the flat behavior of the S — M relation is not 
allowed. This result is of great importance for two reasons 
that we will describe in the following (see section 3.2). 

In their paper, CT10 conclude that results for their ref- 
erence model argue against the universality of the column 
density proposed in D09. However, the same authors showed 
that describing the DM halo adopting a Burkert model and 
a Chabrier IMF (similar to those used in D09), one obtains 
that 

Lv 



log<S(1.6r s ) = 0.02 log( 



10 n i G 



+ 2.65 



(11) 



implying that 5(1.6r_g) is constant with Ly but using a 
Salpeter IMF they got S oc Ly 20 , arguying against univer- 
sality and for a strong model dependency. 

In Fig. 2, we plot S versus A/200 for masses in the range 
10 10 -10 14 ' 5 A/©. The plot uses the same symbols of Fig. 1, 
but now for clarity we did not plot the SIM model of B10, 
the best fit linear relation using the D' Agostini (2005) fit. 
The plot shows that CT10 fit gives larger values than the 
Maccio et al. (2008) and B09 N-body simulations. The S-M 
relation obtained from the present paper shows a flattening 
for masses smaller than 5 x 10 10 Mq. This flattening is the 
result of the fact that for dwarf-galaxies the density profiles 
tend to become shallower and fit with Burkert 's profile are 
better than those obtained through NFW model. This agrees 
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Figure 3. a poH) in units of Mq/pc 2 as a function of galaxy mag- 
nitude for different galaxies and Hubble Types. Data are from: 
Spano et al. (2008) data (empty small red circles), the dwarf 
irregulars (full green circles) N 3741 (M B = -13.1) and DDO 
47 (Mb = —14.6), Spirals and Ellipticals investigated by weak 
lensing (black squares), dSphs (pink triangles), nearby spirals in 
THINGS (small blue triangles), and early-type spirals (full red tri- 
angles). The long dashed line is the result of D09. The solid line 
is the result of the present paper, b < S >o^bary an< i 9bary(ro) 
are plotted as a function of the B-band absolute magnitude of 
the galaxies. From the original sample of D09, here, as in G09, 
are used the dwarf spheroidals data, NGC 3741, DDO 47, and 
the two samples of spiral galaxies (Spano et al. 2008; de Blok 
et al. 2008) which all together span the whole magnitude range 
probed in the original sample (D09). The magenta triangles are 
the dwarf spheroidal galaxiesl, the left green point is NGC 3741 
and the right green point is DDO 47, and the empty red circles 
and filled blue triangles symbols are the Spano and THINGS (de 
Blok et al. 2008) spiral galaxies samples, respectively, c Same as 
panel a, but in this case our po r O is calculated without taking 
account of baryons in our SIM, as in dissipationlcss simulations. 



with the claim of D09 and G09. However, we want to stress 
that in the case of D09 and G09, because of the large error 
bars the claim of constancy of the surface density for all 
masses is not univocal. It is very easy to show that the points 
in D09 and G09 can be easily fitted by a line with slope 
of 0.16 (5" oc A/ 016 ). Comparing the results with N-body 
simulations (B09 and Maccio et al. 2009) shows that for 
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larger masses in the cluster mass region, the prediction for 
S of the present paper are in agreement with simulations but 
for smaller masses tend to diverge. The reason of this is due 
to the fact that N-body simulations uses just CDM particles 
and do not take account of the effects of baryons and the 
others taken into account in the present paper model. 

In order to compare our results to those of D09, and 
G09, we need to find a way to express our result in terms of 
Pod = poro- 

In the following, we show some examples for known pro- 
files. Before we write the S(r) for them. For the pseudo- 
isothermal profile one obtains: 



Siso(R) = 



2np c r c 
R 2 



VR 2 + rl - 



For the NFW density distribution, we have 



Snfw(R) = 



R 2 



arctan -J R 2 /r 2 — 1 

y/R?M - 1 



+ 



(12) 



.(13) 



Analytical relation between parameters of several profiles 
can be obtained fitting the same rotation curve (B09, CT10). 
For an ISO and NFW profile, one can take an ISO profile, 
calculates the relative rotation curve and fits the result with 
a NFW profile. One find: r 3 ~ 6.1r c , p s — 0.11r c . Compar- 
ing the column densities for NFW and ISO profiles, whose 
parameters are related through the relations now written, 
one obtains: 

SNFw{r B ) 



>(6r c ) 



i 0.91 



Similarly for the Burkert profile and NFW: r s 
p s ~ 0.37/9Q and 



SNFw(r s 



0.! 



(14) 
1.6ro, 



(15) 



Ss«rfcert(X.66ro) 

The difference between the column densities Snfw, 
S Burkert and Si so turns out to be less than 10% (CT10). 

Then S N Fw(r s ) ^ 0.98S B urkert{l-&r ) « 1.89r po- 

In Fig. 3a, Fig. 3b, and Fig. 3c, we compare the result 
of the present model with D09 and G09 results. Fig. 3a plots 
poro in units of Mo/pc 2 as a function of galaxy magnitude 
for different galaxies and Hubble TypeQ The quoted fig- 
ure shows that the column density is indeed constant for 
masses smaller than ~ 5 x 10 10 Mq and luminosities smaller 
than Mb ~ —14. However, a systematic change of the aver- 
age column density poro as a function of the object mass is 
clearly present for larger masses, the data being well fitted 
as poro tx M a with a = 0.16 ± 0.05. Such results make, 
therefore, us safely argue against the universality of Sdm 
(i.e., Sdm = const) claimed by D09. 

Fig. 3b shows instead the stellar column density de- 
noted here (E)o& to be consistent with GO£0- For this same 

4 The empty small red circles are taken from Spano et. al. (2008). 
The URC (solid blue line) , the dwarf irregulars (full green circles) 
N 3741 (M B = -13.1) and DDO 47 (M B = -14.6), Spirals and 
Ellipticals investigated by weak lensing (black squares), dSphs 
(pink triangles), nearby spirals in THINGS (small blue triangles), 
and early- type spirals (full red triangles). The long dashed line is 
the result of D09. The solid line is the result of the present paper. 

5 Prom the original sample of D09, here, as in G09, are used 
the dwarf spheroidals data, NGC 3741, DDO 47, and the two 



reason, we also show the corresponding Newtonian accel- 
eration, gb{ro) = 7rG(E)o6 in cm/s 2 . Both quantities are 
plotted as a function of the B band absolute magnitude us- 
ing the same dataset in G09 which is obtained from the 
D09 sample removing some problematic objects (see G09 
for details). As for Sdm, we again find that (E)ot. is con- 
stant for M < 5 x 10 10 M , while it scales with M as M a 
with a = 0.18 ± 0.05 for larger masses. We therefore agree 
with both D09 and G09 at low masses, but the slope at high 
M is a clear evidence against their claim that the baryons 
acceleration at ro is constant over the full mass range. 

In Fig. 3c, we plot the surface density of dark matter for 
a system composed only of dark matter, as in dissipationless 
simulations. In order to get the surface density just of dark 
matter, we did not take into account baryons in our SIM. In 
this case, we get p r oc M a20±a05 . 

As reported in the introduction, also NRT10 questions 
D09 and G09 results, who claimed a universal DM projected 
density for all galaxy types and masses (see also Kormendy 
& Freeman 2004). They reproduce their results reasonably 
well, finding a fairly constant surface density on average 
across a large mass range of dwarfs and late types, but also 
suggesting a systematic increase with mass as found by B09. 

One objection that could arise at this point is that the 
results presented in Fig. 3, are not better fits to D09, and 
G09 data than a costant value. Our result is not a fit to 
the data, but it is the result of a theoretical model, which 
predicts, for example in the case of the surface density of 
baryons a = 0.18 ± 0.05, and equal error bars in the other 
cases. Our result with an error of ±0.05 and not ±0.2 (as in 
D09, and G09), differs at a level of 3 a from G09 (and D09) 
result, and specify among all the slopes permitted by G09 
([—0.2,0.2]) a smaller set. This has a strong physical conse- 
quence, namely that the surface density is not a universal 
quantity. In fact, in order the acceleration scales gi(ro) (or 
the surface density), (with i = {DM, bar}) was an univer- 
sal quantity as argued by D09 and G09, one should find no 
correlation with any stellar or DM property. The wide scat- 
ter in the slope obtained by the authors (0 ± 0.2) prejudice 
this conclusion, since it implies at least a correlation with 
Mb- D09, and G09 results are compatible with a flat slope 
implying no correlation, but are also compatible with a non 
flat slope implying correlation with Mb- Our result is not 
compatible with a flat slope. 

Nevertheless the previous discussion, one could argue 
that even if the model is consistent, it is not necessarily 
more correct than the constant. 

In order to further clarify the previous issue, we re- 
peated G09, and D09 calculations in a similar way they did, 
using a Burkert's profile fit, but for a larger sample. The 
details of the analysis: the data used, the way the fits were 
performed, were published in Cardone & Del Popolo (2012). 
Here, we summarize the main features of the paper. In the 
quoted paper, we searched the literature for systems with 



samples of spiral galaxies (Spano et al. 2008; de Blok et al. 2008) 
which all together span the whole magnitude range probed in 
the original sample (D09). The magenta triangles are the dwarf 
spheroidal galaxiesl, the left green point is NGC 3741 and the 
right green point is DDO 47, and the empty red circles and filled 
blue triangles symbols are the Spano et al.(2008) and THINGS 
(de Blok et al 2008) spiral galaxies samples, respectively. 
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high quality rotation curve data probing the gravitational 
potential with good sampling, large radial extension (i.e., 
up to R > R pt = 3.2Rd with Rd the disk scale length) and 
small errors. In order to reduce the potential impact of beam 
smearing on the circular velocity, we used, whenever possi- 
ble, Ha data or a combination of HI and Ha. Our sample 
contains 58 galaxies (also containing G09 and D09 galaxies) 
from different sources (de Blok & Bosma 2002; Simon et al. 
2005; Spano et al. 2008; THINGS; Oh et al. (2008); Swa- 
ters et al. 2011) spanning nine orders of magnitude in the B 
band luminosity (roughly, —22 ^ Mb ^ —13) and different 
spiral classes (from dwarfs system to low and high surface 
brightness galaxies). In order to be consistent with G09, we 
assumed a Burkert model to fit the rotation curve. 

The model parameters were obtained through a likeli- 
hood analysis by maximizing £(p) oc exp [— x 2 (p)/2] with 



2 

X = 



Af r 

E 



obs i \ 

v c (n) 



l (n, P ) 



(16) 



where v° (ri) is the measured circular velocity for the 
i-th point (with at the error), w* (^iiP) the theoreti- 
cally predicted value for the given set of parameters p = 
(T*, V v ir, Cuir) and the sum is over the M data points. A 
Markov Chain Monte Carlo (MCMC) method running three 
chains checking the convergence through the Gelman & Ru- 
bin (1992) criterion was used to efficiently explore the quoted 
3D space of parameters. 

After having constrained the halo model parameters 
and the disc stellar M/L ratio, the Newtonian acceleration 
of both the baryon (disc + gas) and DM components can 
be estimated through Eqs. (1) and (4) in Cardone & Del 
Popolo (2012). The uncertainties on the fitting parameters 
and the correlations among them was obtained evaluating 
gi{rj) along the Markov chain for each galaxy and use the 
resulting values to estimate the median and 68% confidence 
ranges. Note that, here, the label i denotes the galaxy com- 
ponent (baryons or DM), while j refers to the reference ra- 
dius adopted. 



Our method has two differences when compared to G09. 
Firstly, we fit the full rotation curve with the Burkert model 
thus being also able to check whether or not the model well 
reproduces the observed data. G09 did actually not perform 
any fit to the rotation curve data, but used literature data 
and an approximate law to convert the scale length of the 
model used by other authors to their Burkert model ro ra- 
dius. This choice does not automatically ensure that the 
theoretically predicted rotation curves well match also in 
the outer regions. However, the scaling law used by G09 in- 
deed provides values of gi(ro) consistent with our ones with 
(Alog5DAi(r )} = -0.03 and (A log g bar (r )) = -0.02 as 
estimated by the 18 galaxies common to both sample^ 

A second difference with G09 concerns the estimate of 
the uncertainties. As previously reported, we have deter- 
mined (rather than assumed) the halo model parameters 
and the stellar M/L ratio, and this allowed us to attach 
to each ffi(ro) value an error fully taking into account both 
the fitting parameters uncertainties and their correlations. 
We then add in quadrature a systematic error related to the 
uncertainties on the method used to determine gi(ro). This 
was obtained comparing with the G09 values for the galaxies 
common to both samples and add in quadrature 0.23 to the 
statistical errors on loggi(ro) (see Cardone & Del Popolo 
(2012) for details). 

In order to check the universality hypothesis, we fit the 
following log -linear relations 



log g t (r ) = a log 



10 11 M Q 



+ b 



log#j(ro) = aMs + b 



using the Bayesian method, getting (we show here just the 
gtar - M B case) for g bar - M B , a = O^Sl+g 1 ;^, b = 
13.04lo;°3g, and the intrinsic scatter a int = 0.23+°;^ ob- 
tained through the DAgostini (2005) Bayesian method. 



We miss 10 galaxies from the G09 sample because of different 
reasons. 
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In Fig. 4, we plot our data (black) superimposed to G09 
(red) and the best fit line. 

Fig. 4 shows that in the range —20 ^ Mb ^ —15, 
our data are well superimposed to the G09 ones except 
for few cases where our estimates is larger than that in 
GOSfl, and that if we exclude the galaxies with log gb ar (ro) > 
-8.0 and add the Milky Way dwarfs with Mb > -12, a 
constant gbar{ro) would be preferred. However, the situa- 
tions changes when taking into account the four systems 
with logg bar (r ) > -8.0 (namely, NGC2841, NGC3521, 
NGC4736, NGC6946), which belong to the THINGS sam- 
ple so that they have a high quality rotation curve and a not 
small inclination (as estimated from the HI data) . These sys- 
tems were not considered in G09. The reason for the discrep- 
ancy on the constancy of gb ar (ro) with Mb is mainly due to 
our sample including galaxies with larger acceleration val- 
ues. The previous result clearly shows the non-universality 
of the quoted relation. 



4 CONCLUSIONS 

In the framework of the SIM model of DP09, in this pa- 
per the correlation between the central surface density and 
the halo core radius of galaxy, and cluster of galaxies, dark 
matter haloes was analyzed. The main result is that r* 
(the halo characteristic radius)is not an universal quantity 
as claimed by D09 and G09. We obtained a surface den- 
sity , S oc M 016±0 - 05 , and S oc M°- 18±0 05 , for masses 
> 5 x 1O 1O M , for the cases examined in D09 and G09, 
respectively. Our results suggest that the constant surface 
density, claimed by the previously quoted authors, is prob- 
ably unlikely. 
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